arXiv:1507.00080vl [math-ph] 1 Jul 2015 


On the Attractor for the Semi-Dissipative Boussinesq Equations 

Animikh Biswas'^’*, Ciprian Foias’’, Adam Larios*^ 

°‘Department of Mathematics & Statistics, University of Maryland, Baltimore County, 1000 Hilltop Circle. 

Baltimore, MD 21250, USA 

^Department of Mathematics, Texas A&M University, College Station, TX 77843, USA 
‘^Department of Mathematics, University of Nebraska-Lincoln, Lincoln, NE 68588-0130, USA 


Abstract 

In this article, we study the long time behavior of solutions of a variant of the Boussinesq 
system in which the equation for the velocity is parabolic while the equation for the tem¬ 
perature is hyperbolic. We prove that the system has a global attractor which retains some 
of the properties of the global attractors for the 2D and 3D Navier-Stokes equations. More¬ 
over, this attractor contains inhnitely many invariant manifolds in which several universal 
properties of the Batchelor, Kraichnan, Leith theory of turbulence are potentially present. 
Resume. Dans cet article nous etudions le comportment an temps inhni des solutions d’une 
partiellement dissipative system du Boussinesq, dont une est parabolique et I’autre est hyper- 
bolique. Dans ce but nous introduire un attracteur universel qui retient plusieurs proprietes 
des attracteurs universels des equations de Navier-Stokes en dimension deux on trois, qui, in 
particulier, contient une inhnite de varietes invariantes dans lesquelles plusieurs proprietes 
universelles de la theorie de la turbulence bidimensionnelle de Batchelor, Kraichnan et Leith, 
sont potentiellement presentes. 
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1. Introduction 


A critical part of understanding turbulence in geophysical flows lies in understanding the 
large-time dynamics of the system. A tool that provides deep insight in studying large-time 
dynamics of certain systems is the study of global attractors, which typically requires the 
underlying dynamical system to be dissipative; see for instance [3,0, 39, and references 
therein. For instance, in case of the fully dissipative Boussinesq system, the study of the 
long term dynamics via global attractors has been accomplished in |20(]. However, in certain 
physical regimes, the dynamical system governing geophysical flows is modeled to be only 
semi-dissipative. That is, the system is dissipative in some variables, but not in others. This 
means that the usual notion of global attractors no longer applies. In this work, we propose 
a new type of global attractor—which we call a weak sigma-attractor—that captures the 
large-time dynamics of a well-known semi-dissipative model for ocean flows known as the 


semi-dissipative 2D Boussinesq system [3^. The notion of the weak sigma-attractor, which 


in some respects is a hybrid of the attractors of the 2D and 3D Navier-Stokes equations 


25|, exploits the fact that while only part of the system is dissipative, the other part 


has a hyperbolic structure that gives rise to certain conserved quantities which in turn 
correspond to invariant sets in the phase space. We expect that this approach, which is new 
to the best of our knowledge, will be useful in studying the long-term dynamics of many 
other semi-dissipative systems. 

The global attractor for a dissipative dynamical system is usually deflned as either the 
maximal compact invariant set, the minimal set which uniformly attracts all bounded sets, or 
the set of points on complete bounded trajectories. For many dissipative systems (e.g., the 
heat equation, the 2D Navier-Stokes equations, the ID Kuramoto-Sivashinsky equations, 
and many others), these deflnitions are equivalent. For the semi-dissipative Boussinesq 
system, this equivalence fails, and therefore the deflnition must be modifled. In particular. 
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compactness must be sacrificed. However, as we will show, it can be recovered in a certain 
local sense, as the attractor we propose will be shown to be a countable union of weakly- 
compact invariant sets that retain several properties of the attractor for the 2D Navier-Stokes 
equations. Moreover, the attraction properties we describe will be in the sense of the weak 
topology rather than the str ong topology. This aspect is similar to the situation for the 3D 
Navier-Stokes equations (23|, [2^, where, due to the lack of knowledge about the global well- 
posedness of the system, one studies a weak attractor rather than a strong global attractor. 
However, the system we study is known to be globally well-posed, and thus, the use of the 
weak topology is intrinsic to the semi-dissipative nature of the system under consideration. 
In contrast to this, if the global well-posedness for the 3D Navier-Stokes equations were to 
be established, the corresponding weak attractor will immediately become an attractor in 
the usual (strong) sense. 

A curious feature of the 2D semi-dissipative Boussinesq system (HU pertaining to the 
study of 2D turbulence deserves mention here. In particular, the equations potentially 
combine a wealth of turbulent dynamics in a single system. This is due to the fact that 
the right-hand side of the momentum equation is a time-dependent dynamical variable, 
but that its norms are time-independent (1 < p < cxo). Thus, the velocity component 
of fll.ljl will indeed contain turbulent dynamics corresponding to all Grashof numbers, so 
long as one can establish that arbitrarily large effective Grashof number-defined as the 
magnitude of the divergence-free part of the force that drives the velocity equation-can 
be achieved for large time. We discuss this issue in greater detail in Section [HI On the 
invariant subsets mentioned above, the system is expected to have time statistics consistent 


with the Batchelor-Kraichnan-Leith |J, |35| empirical theory of 2D turbulence, albeit with 
different Grashof numbers varying from subset to subset. We can therefore think of the 2D 
semi-dissipative Boussinesq system as a platform which hosts an enormous variety of 2D 
turbulent dynamics. It appears to be a model system for the study of attractors, in the 
same way that G. Elegans and Drosophilia are model species for study in biology. (Similar 
statements may be made in the 3D case.) We emphasize that this system was not created 
in an ad-hoc fashion to have these particular properties, but instead arises naturally as a 
relevant system in geophysics. We will show that for system fll.ip . the weak sigma-attractor 
is not only non-trivial, but it is an extremely rich proper subset of the phase space. We 
hope that the study of the long-term dynamics of this system will inspire further research 
leading to greater insights into 2D turbulence, and perhaps suggest new approaches to this 


area. 


The paper is organized as follows. In Section [LTl we describe the semi-dissipative Boussi¬ 
nesq system, which is the main subject of our study. In Section m we lay out notation and 
preliminary material. In Section 12.11 we recall well-known well-posedness results which will 
be useful in the sequel. In Section |3l we prove that smooth solutions to system fll.ip have 
the backward-uniqueness property. In Section jH we propose the notion of a weak sigma- 
attractor for system da, describe several new classes of steady states (as well as time 
variant solutions) of the system, and state the main theorem on the properties of the weak 
sigma-attractor. In Section |H1 we prove the properties discussed in Section |U In Section [HI 
we discuss connections to turbulence as evidenced in the 2D Navier-Stokes equations with 
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forcing in all Fourier modes, while in Section [TJ we discuss the projection of the attractor 
onto the energy-enstrophy plane. We hnd certain points in this setting which correspond 
to steady columnar flows which are in hydrostatic balance. In Section IH we discuss some 
open questions related to this study. Section [9] is an appendix, where some standard facts 
are collected for reference of the reader. 

1.1. Semi-dissipative Boussinesq system 

The semi-dissipative 2D Boussinesq system (without rotation) in the periodic domain 
D := [0, for time t > 0 is given by 


dju -I- (u ■ V)u -I- Vp = uAu -I- dg. 

in D X R_|_, 

(l.la) 

dtO + (u ■ V)d = 0, 

in D X M_|_, 

(1.1b) 

V ■ u = 0, 

in D X R+, 

(1.1c) 

u(x, 0) = Uo(x), d(x, 0) = do(x). 

in D, 

(l.ld) 


equipped with periodic boundary conditions in space js^. Here z/ > 0 is the fluid viscosity. 
The spatial variable is denoted x = {xi,X 2 ) G D, and the unknowns are the fluid velocity 
held u = u(x, t) = (mi(x, t), M 2 (x, t)), the huid pressure p(x, t), and the function 9 = 6'(x, t), 
which may be interpreted physically, e.g., as the temperature variable. We write g = (0, g)"’" 
for the constant, upward-pointing gravity vector, where g is the (scalar) acceleration due to 
gravity. It is straight-forward to show the Galilean invariance and mean-preserving property 
of fll.ip (see Section IT^ . Therefore, we henceforth assume that f^udx = 0 and f^9 dx = 0, 
and similarly for the initial data. Furthermore, all function spaces used in this work will 
be assumed to contain only mean-free elements, unless otherwise indicated. We remark 
that, with minimal adjustments, the results here are also valid in the presence of a Coriolis 
rotational term. 


2. Notation and Some Specific Preliminaries 

Let J- be the set of all vector-valued trigonometric polynomials with periodic domain 
D = := M^/(LZ^) = [0, L]^. We dehne a space of test functions 


V:=<(pGJ^;V-(p = 0 and / (p(a;) dx = 0 > . 


/T2 

We denote by and El™' the usual Lebesgue and Sobolev spaces over T^, and dehne H and 
V to be the closures of V in and EI^ respectively. We will often use the notation (•, •) to 
denote pairs, so to avoid confusion, we denote the inner products on H and V respectively 
by 

duj dvi 


(u,v):=V/ UiVidx and ((u,v)):=V ^ ^ 

and the associated norms, respectively by 

||u||i2 := (u, , ||u||hi := ((u, . 


dx. 
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This notation will also be extended to apply to tensors, such as Vu, in the natural way, which 
should not be a source of confusion. Note that || • || is a norm due to the Poincare inequality, 
dll. below. We denote by V the dual space of V. The action of V' on V is denoted by 
(•, ■)y, y. Note that we have the continuous embeddin gs V ^ H ^ V. Moreover, by the 
Rellich-Kondrachov compactness theorem (see, e.g., [l|, ll5[), these embeddings are compact 
due to the boundedness of the domain T^. 

We denote hy ^ H the Leray-Helmholtz projection operator (i.e., the orthogonal 

projection onto divergence-free vector spaces), and the Stokes operator A := —Po-A with 
domain 'P(A) := fl R. In our case of periodic boundary conditions, it is well-known 


that A = —A (see, e.g., [8|, |^). A ^ : P —)■ P is a positive-dehnite, self-adjoint, compact 


operator, and therefore has an orthonormal basis of eigenfunctions (pk corresponding to a 
non-increasing sequence of eigenvalues (see, e.g., @,11^). We label the eigenvalues Xk of A so 
that 0 < Ai < A 2 < As < • • ■. Furthermore, for all w G Id, we have the Poincare inequality. 


«;o||w||l2 < ||w||hi, kq\=2'k/L. 

Due to fl2.ll) . for w G P(A), we have the norm equivalences 

\\Aw\\l 2 = ||w||h 2 and ||Vw||i 2 = ||w||hi. 
In fact, for any s G M, one can show that 

||w||h= ^ ||.4“^^w||i2 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


for w G where P(y4®A) and are dehned in a natural way by the eigenvalues 

and eigenfunctions of A (see, e.g., [sl, 40 


We also note the Ladyzhenskaya inequality, which in 2D reads 




s 

(2.4) 

|u|Il2^I|Au||^(^ 

(2.5) 

P<,((wi ■ V)W2) 

(2.6) 


We use the standard notation 


for wijWoG V. We note several important properties of B which are proven, e.g., in 


19l. l40l . l42l | . The operator B dehned in fl2.6p is a bilinear form which can be extended as 


a continuous map B : V x V ^ V. Furthermore, for Wi, W 2 , W 3 G V, 

(P(wi, W2), Ws)^,^^ = - (P(wi, Ws), W2}y,y , 
(B(Wl,W2),W2}y,y = 0. 
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(2.7a) 

(2.7b) 















The following inequalities hold for smooth functions: 


I (5(wi, W2), W3) I < C'||wi||^'^2^||wi||J/i^||w2||^'^2^||w2||J/i^y^ (2.8a) 

I (5(wi,W 2), W3) I < C'||wi||^'^2^||wi||J/i^||w2||Hi||w3||^i^||w3|li^^^ ( 2 . 8 b) 

I (5(wi, W2), W3) I < (T||wi||^'^2^||Awi||^(^||w2||hi||w3||2,2 (2.8c) 

I (B(wi, W2), W3) I < C'||wi||2.2||w2|l][/i^||A^/V2||^'^2^||w3||i2^ (2.8d) 

These results also hold if B is extended to allow W2 and W3 to be scalar-valued functions in 
the indicated spaces (we still require V ■ wi = 0), and A is replaced by —A. 


2.1. Known Well-Posedness Results 

^he global well-posedness of system fll.ip in the case (uq, Oq) G (H^nld) x was proven 

The requirements on the initial data 
see also jl^, [s^ for related 


m 


27 


in j^, and in the case (uq, 6*0) G (H^ fl Id) x i 
were weakened in j^, where the following theorem was proven 
results and additional discussion). 


Definition 2.1. We say that a pair (n, 6) solves fll.ip if (n(f), 6{t)) E H x for each f > 0, 
u and 6 are both periodic and mean-free, and they satisfy 


{u{s), {s)) ds + / (n(s) (8 n(s), V<h) ds 


(2.9a) 


= (uo,$(-,0)) - z/ / ((u(s),$(s))) ds+ / (d(s)g,<l)(s)) ds, 


{9{s), (j)'{s)) ds+ {u{s)e{s), V0(s)) ds = (do, 0(-, 0)) 


( 2 . 9 b) 


Jo Jo 

for all mean-free, space-periodic scalar test functions (j){x,t) G x [0,T]), such that 

(f){x,T) = 0; and for all mean-free, space-periodic vector-valued test functions ^{x,t) G 
[C'°°(fl X [0,T])]2 such that V ■ <!>(•, f) = 0, <!)(•, T) = 0. 

A pair (u(f), 0(t)) E H x f G M is a global solution if 


(n(s), <h'(s)) ds-|- / (u(s) ® u(s), V<h) ds 

J —00 

>D noo 

= -V I ((n(s),<l)(s))) ds / (d(s)g,<l>(s)) ds, 

X) J —00 

POO 

(d(s), 0'(s)) ds-F / (n(s)d(s), V0(s)) ds = 0, 


' —00 

POO 


J —00 
'*00 


(2.10a) 


(2.10b) 


where 0 G C'°°(r2 x (— oo,cxd)) and <h G [C'°°(fl x (— cx), cxd))]^ are exactly as before, except 
that they are assumed to be compactly supported in time. As mentioned before, the pressure 
p is recovered via fl2.20p . 
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Theorem 2.2 ([^). Suppose (uqj^o) G H x Then, for any T > 0, the systen^ fll.ip has 
a unique solution (u, 6) such that 


ueC{[0,T],H)nL^{[0,T],V), 

^eLH[0,T],V% 

If additionally, uq G V, then 

ueC{[0,T],V)nL\[0,T],V{A)), 

-eL\[0,T],H), 


9eC{[0,T],L^), 

^GL^([0,T],H-3/2). 





Recently, results on well-posedness and persistence of regularity were improved including 


the treatment of the bounded domain setting [28|, l29j. The following two theorems collects 


some of the results of the aforementioned works. 


Theorem 2.3 (Combined results from [1,0, H, 28, 2^). Let s > 0. Suppose Uq G 
and 00 G ff. Then u G ^([0, cx)), ei+^) n Ll^{[0, oc),B^+^) and 9 G ^([0, oo), H^). 


Theorem 2.4 (Combined results from 0,0, m, [H, H^). Let s > 1. Suppose uq G fl V, 
and 00 G Then u G ^([0, cxd), H^) fl L|,j,([0, cxd), and 0 G ^([0, cxd), H^). 


Next, we recall the well-known result (see, e.g., H) from which it follows that the 
transport equation has unique solution which preserves all L^ norms as well as distribution 
functions. Recall that for a function ip G L‘^{fl), the corresponding distribution function 
: R —)■ [0,1] dehned by 




1 

M 



X(-oo,p](</2(a;)) dx 


1 

M 


|{a; G : (p{x) < p}| . 


Proposition 2.5. Let v G L^(0,T; V) and 9 q G Lf{VL), the Cauchy problem 


( 2 . 11 ) 


ra,0 + V-(v0) = O, 

(2.12a) 

\ 0(x,O) = 0o(x), 

(2.12b) 

has a unique solution 9 = 9{'x,t) such that for a.e. t >0 and all p G R, 


Fe{-,t){p) = Feoip). 

(2.13) 

In particular, for all p G [1, cxo], and a.e. t > 0, 


1 0(^)1 LP = 1 0o| LP- 

(2.14) 


^in the appropriate weak formulation, (EH) below 


7 














The proof is well-known, but we sketch it here for completeness. 


Proof. The existence and uniqueness is well-known; see, e.g., and the references therein. 
Equation fl2.13p is also well-known, and follows by noting that, since v G L^(0,T; V) and 9 
satisfies fl2.12ap . 6 is constant along streamlines. A change of variables (sometimes called the 
“back to labels” map), using the fact the V ■ v = 0, then yields fl2.13p (see, e.g., 0, [f^ for 
more details). Furthermore, fl2.14p holds for p G [1, cxo) by applying fl2.13p to the well-known 
identity (see, e.g., j^, VI. 10) 


p 

LP 


ffdF^ip). 


(2.15) 


The case p = oo is also classical, and can be proven, e.g., by using the Hopf-Stampaccia 


technique (see, e.g., |30|]). Namely, one notices that Q := 6 — ||6 'o||l°o also satishes fl2.12ap . 
and then takes the inner product of the resulting equation with 0+ and 0“, the positive 
and negative parts of 0. Integrating by parts and using the divergence-free condition yields 
||0||i = O. □ 


2.2. A Functional Form of the equation 

As is customary for the Navier-Stokes equations, applying P„ to fll.lap (see Section [2] for 
notation) yields the functional equation 


— + uAu + B{u, u) = Pa{0g). 

It is straight-forward to check that 

p,{eg) = »g - 9V(4 ,a-‘») = 

As this is an orthogonal decomposition, 

||0g||i. = ||P.(0g)||i.+/||V(4,A-^0)||i.. 


Thus, 


g-^\\PAm?L- = ll«ni» - (VA-‘S„«, VA-‘S„«) 
= ||«||i.+ (AA-‘a„».A-‘9„«) 

= wnh - iK-A)-‘'=9„«iii, 

= ll»lli. - \\RA\h = \\RA\h. 


(2.16) 


(2.17) 


(2.18) 


(2.19) 


where Ri := (-Aj-'/’clj,,•! = 1,2 are the Riesz-transforms, which can also be defined by 
their Fourier symbols Pj(k) := i/cj/|k|, z = 1, 2 at wave number k = (/ci, k 2 ). 




















Due to the mean-free assumption, if one obtains a pair (u, 6) which solves 02.161) and 
Ol.lbD - Ol.ldD . then one can obtain the pressure p by taking divergence in Ol.lal) and noting 
that 


p = - ^ Rij{uiUj) + gA Ry = A 


( 2 . 20 ) 




where is taken with respect to the periodic boundary conditions and the mean-free 
condition. Thus, we will also refer to a pair {u(t),9(t)) solving 02.16P and Ol.lbp - Ol.ldp as 
a solution to the system (Hj. 


3. Backward Uniqueness 

Our main result in this section is that the system 01 .ip has the property of backward 
uniqueness for sufficiently smooth solutions; namely, that if two (sufficiently smooth) solu¬ 
tions agree at some point in time, then they agree for all previous times. While this is known 
for the 2D Navier-Stokes equations and for the transport equation, the result for the coupled 
system 01 .ip does not seem to follow in a straightforward manner from the corresponding 
results for the uncoupled system. Our technique is inspired by that of [^. The backward 
uniqueness result will be useful in studying the asymptotic dynamics of 01.11) . For instance, 
together with well-posedness, it implies that if any two global trajectories intersect, then 
they must be identical. 

Let us recall that the method in (which applies to a wide class of fully dissipative 

equations, including the 2D Navier-Stokes equations) is to consider two solutions ui and U 2 

which are not identical at some time, and then to show that the log of the norm of their 

difference, u := ui — U 2 , remains bounded for all future times, so that the solutions can 

never coincide in the future (this is what is meant by backward uniqueness). The estimates 

in typically require one to show that the Dirichlet quotient X{t) := ||u ll^i/l|5|| I 2 grows 

at most exponentially in time. The lack of diffusion in fll.lbp poses a hurdle in showing 

such a bound in our case. However, we show below that by considering instead the hybrid 
|| 2||2 

quantity I{t) := ygip +g°]|g ||2 > diffusion in fll.lap alone is enough to gain the necessary 
control over both variables u and 9. Note that I{t) is not a Dirichlet quotient. 

Theorem 3.1. Fori = 1, 2, let (uj, 9i) be solutions to fll.ip . such that (uj, 9i) G L°°((Ti, T 2 ); Hf] 
jj3) ^ L°°((Ti, T 2 ); H^). If there exists a time T* G (Ti,T 2 ) such that (ui(T*), 6 'i(T*)) = 
(u2(T*), 6'2(T*)), then for a.e. t G (Ti,T 2 ), we have (ui(f), ^(f)) = (u2(t), 6 * 2 (t)). 

Proof. Denote u := Ui — U 2 , 9 := 9i — 6 * 2 , p := pi — P 2 - 

dtu + (ui ■ V)u (u ■ V)u 2 -h Vp = uAu + 9g, (3.1a) 

dt9 + (ui ■ V)9 + (G ■ V)02 = 0. (3.1b) 

Suppose that there exists some T* > 0 for which both u(T*) = 0 and 9(T*) = 0, but 
that there also exists some to G (Ti,T*) such that u(fo) ^ 0 or 9(to) ^ 0. Due to the 
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continuity in time given by Theorem 12.41 there is an interval [to, to + 5) on which u ^ 0 or 
9^0. Let us denote by [to, ti) the largest such interval, noting that we must have u(ti) = 0 
and 9{ti) = 0. For the remainder of this section, we work on this interval. Taking the inner 
product of fld.lap with u, and with Au, and of fl3.1b|l with 6, integrating by parts, and using 
fl 2 Tbll yields 


^^I|u||i 2 = -I^||u||^i - ((u ■ V)U 2 , u) + (oz, , 

^^I|u||hi = - ((u ■ V)U2, Au) + (eg, Au) , 

1 d 

2 dt' 

Consider the quantity L dehned on [to,^i) by 


\\e\\l. = -(iu-V)9,,9). 


(3.2a) 

(3.2b) 

(3.2c) 


L{t) := log 


Mh + gmh, 


^\\l^ + g^dWl^ 


(3.3) 


Using fl3.2ap and fl3.2cp . we compute 

-1 


dL 

dt 


^Wh + g^ 


(liMhAdimh 


12 I ^2||fli|2^ \2dt" 2 dt' 


jy\ 


u|Ihi + ((u ■ V)U 2 , u) - (eg, uj +g‘^ ({u- V) 6 ' 2 , 9j 


u 


L 2 


+ g^ 


Il 2 


< u- 


u 


Hi 


Mh+gmh 

:= z/J + // + ^ + III 


((U- V)U2,U) ^ 1 ^ 


g^({u-v)e2,e 


u 


L 2 


+ r 


Il 2 


u 


L 2 




IL2 


(3.4) 


where we have used Young’s inequality. Thus, in order to hnd a bound for dL/dt, we search 
for a bound on the quotients I, II, and III in the right-hand side of 03.41) . First, note that, 
due to 02.8bD and 02.8cl) . 


((u- V)U2,u) < C||U2 ||hi||u||li||u||h. < C||U2||Hl||u|||a + ||u||i,,, 

/((a . V)92.h < Cg^l|A92|k.||S||e.||9|U. < C9‘‘||A92||?,2||9||i2 + ||S||^.. 


(3.5) 

(3.6) 


Thus, 


jj ((U ■ V)U2, U) ^ g||u2|lHl||u^i2 ^ 


u 


Hi 


\m.+gmh 


im^+gmh 


Mh + gmh 


^ C'||u 2 ||hi- h/, (3.7) 
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and 


III : = 


\(u-v)02,0; Cg^Ae^w 


u 


iL^ll^lljz + II^ IihI ^ ,^„2|| A a Il2 

h + g^e\\h \Mh + 9^nh 


< 


^ Cg‘^\\A92\\‘j^2 + I 


(3.8) 


In order to estimate J, we compute dl/dt. First, we derive an inequality which will be 
useful in estimating dl/dt. 


\\Ami, - i\\u\g, = \\Am, 


u 


Hi 


U 


L 2 




|2 

Il 2 


>ll^5||i.-2- 


u 


Hi 


r2||~ll2 


U 


Hi 


U 


L 2 


+ 9^ 


|2 

Il 2 


= \\Au\\l2-2I{Au,n)+P\\u\\l2 
= IIAu-Jull^^. 


u 


L 2 


U 


L 2 




|2 

Il 2 


(3.9) 


Using fl3.2|] . we compute 
d 


Idl 
2 dt 


2 II ^IIhI 


1 _d ||~||2 

2ddl^llHl 


U 


|2 1 A( 
I Hi 2 dt ^ 


U 


1,2 




Ii 2 ) 




|2 

Il 2 


U 


L 2 




|2 

Il 2 


U 


L 2 




lL 2 j 


-z/||245ii2 


L 2 


((u ■ V)u 2 , Au) + ( 9g, Au 


u 


L 2 


+ 9^ 


|2 

Il 2 


U 


Hi 


^l|u||ei + ((u-V)u 2 ,u) - ( 6 ig,u) + 5(2 ( (u ■ V) 6 ' 2 , 61 


u 


L 2 




|2 '(2 
Il 2 J 


-i/IIAull^j - ((u ■ V)U 2 , Au) + 


, Au 


u 


1,2 


+ 9^ 


|2 

11,2 


+ /■ 


i^l|u||ei + ((u-V)u 2 ,u) - ( 6 ig,u) + g^ {{u-V)92,9 


u 


1,2 


+ 9^ 


|2 

11,2 


-i/(||Aiiii 2 


L 2 


-^I|u|Ihi) - ((u ■ V)U 2 , Au - Ju) + ( 6 ig, Au - Ju 


u 


L 2 




|2 

Il 2 


+ /- 


g^({n-V)92,9 


u 


L 2 




|2 

Il 2 


Using this with ([2SaD, (El]) and (jM]), we hnd, 


(3.10) 


!(£_ ^ -z/||Au - Jull^a + ||u||hi||Au 2 ||l 2 ||Au - /u ||^2 + 5 (|| 6 I||l 2 ||Au -/u||i 2 


2 dt 


u 


122+^2 
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|2 

Il 2 


(3.11) 
























Mh + gmih 

I|u|ii2 + c/^l|0|li2 2 

^ ^11 ^^211^2 ill 

z/ ^ z/ ^ 2 ’ 

where we have used Young’s inequality for the last estimate. Applying Gronwall’s inequality 
on [to,t), t G [to,ti) now yields 

I{t) < J(0)e^o^W'^^ + - f (3.12) 

^ Jo 

< J(0)e^o'"*-^W‘^* + 1 f 
^ Jo 

where K{s) := _|_ gll^^WllHS / jg bounded on [to, t). By fl3.7p and fl3.8p . this in 

turn shows that II and III are also bounded, so that fl3.4p implies that dL/dt is bounded 
on [to,t). Integrating in dL/dt in t over [to,T*), we hnd that L is bounded on [to,T*). In 
particular, 

(l|u||i2 + 5'^||0||i2j (3.13) 

is bounded on [to,T*), contradicting the assumptions that u(T*) = 0 and 6{T*) = 0. □ 

The global well-posedness results (Theorem 12.2112.4p imply that fll.ip has a well-defined 
semigroup operator S{t) defined for t > 0 by 

S'(t)(uo,6'o) = (u(t),6'(t)). (3.14) 

The backward uniqueness result implies that the semigroup is injective on at least the smooth 
portion of the set which we will define to be the attractor. This in turn implies that S{t) 
can be extended to hold for negative times on the smooth trajectories in that set. 


4. An Adequate Notion of an Attractor 


We look for a notion of an attractor which applies to fll.ip . While some analogies with 
the attractor, ^ns, of the 2D Navier-Stokes equations can be made, there are also striking 
differences. Therefore, for comparison, let us recall that for the two-dimensional Navier- 
Stokes equations, 

A + „^u + B(u,u) = / (4.1) 


(where / G if is time-independent), the global attractor ^ns 
definitions (see, e.g., iH): 


has the following equivalent 
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(I) 

( II ) 


^Ns is the minimal set which uniformly attracts all bounded sets m. H as t ^ oo. 
^Ns is the maximal compact set of all uq E H such that the solution of fl4.ip satisfying 
u(0 ) = uq exists for all f G M (i.e., backward and forwards in time) and satishes 


sup 

teR 


u 


L2 < OO 


(4.2) 


(III) ^Ns is the maximal bounded set X C H such that S^s{t)X = X for all t G 
S'Ns(t) is the semigroup operator for fl4.ip . 


where 


Moreover, it is well-known that the attractor for the two-dimensional Navier-Stokes equa¬ 
tions has hnite fractal dimension (and therefore hnite Hausdorff dimension) (see, e.g., i0). 

Clearly, any reasonable notion of attractor of a system must include the steady states of 
the system. While the above dehnitions are equivalent ways to dehne the global attractors for 
many dissipative systems, this is not the case for system fll.ip due to the lack of dissipation 
in fll.lbp . Indeed, we will describe below (see Remark 14.ip a set of steady states which are 
neither bounded nor hnite-dimensional. Next we describe special classes of solutions of fll.ip 
which play an important role in our study. 


4-1- Special Classes of Solutions and Steady States 

Here we will provide three families of explicit solutions of fll.ip corresponding to purely 
vertical (V), purely horizontal (H) and plane wave solutions. 


(i) Vertical solutions: Let G L^([0,L]) be an arbitrary periodic, mean-zero 

function depending only on the vertical variable: 9^ = 9^{x2). Set = 0, and dehne 
uf = uf{x 2 ,t) by the following unforced diffusion problem: 

IdtuY - = 0 , 

[una; 2 , 0 ) =a^{x2), 

along with periodic boundary conditions on [0, L] and the mean-free condition. Let 
be dehned up to an arbitrary constant bjo 

dx2P^ = ■ ( 4 . 4 ) 

It is easy to check that {ui,U2,0,p) = {{uf,U2),0^,p^) is a solution to fll.lap - fll.lcp 
with initial data mi(x, 0) = a^{x 2 ), M 2 (x, 0) = 0, and ^(x, 0) = 9^. 

Steady State: It is easy to see that the vertical solution (u'^, 9^) dehned above 
converges to the steady state (0, 9^) as t —)■ oo. 

(ii) Horizontal Solutions: Let a^,9^ G L^([0,L]) be arbitrary periodic, mean-zero 
functions depending only on the horizontal variable: 9^ = 9^{xi). Set = 0, 


^Property (14.411 is often referred to as hydrostatic balance in the geophysical literature. 
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= 0, and let be the unique solution of the following forced diffusion 

problem: 

I dtu^ - = ge^, 

=a^{xi), 

with periodic boundary condition on [0, L] and the mean-free condition. It is easy to 
check that {u^, 6 ^:= {{u^,U 2 ), 6 ^,p^) is a solution to fll.lap - fll.lcp with initial 
data 


Mi(x, 0) = 0,M2 (x, 0) = a^{xi) and 6*(x, 0) = 9^. (4.6) 

Since these flows are independent of X 2 , they can be thought of in a geophysical context 
as columnular flows. In Section (jT]), we show that the steady states corresponding to 
these flows also arise as solutions lying on the boundary of a certain set containing the 
attractor after projecting into the Energy-Enstrophy plane. 

Steady State: Let U 2 = U 2 {xi) be the unique, mean-free periodic solution of the 
equation 


= -gd^{xi). (4.7) 

One can check that the horizontal solution (u^, 9^) dehned above converges, as t —)■ oo, 
to the steady state 


Ml = 0, M2, 9 = 9^. (4.8) 

(iii) Plane Wave solutions: We also consider certain plane-wave (PW) solutions, hrst 
introduced in the context of the Navier-Stokes equations in [^. Specihcally, let / = 
f{zfl) and h = h{z) be smooth mean-free function^ which are L—periodic in z, and 
§x a non-zero vector k = (fci, ^ 2 ) £ where ki + k 2 = 0. Below, in fl4.13p . we will 
require / to satisfy a particular heat equation. Dehne u and 9 by 

uP^(x, t) := (/(k ■ X, t), /(k ■ X, t)), (4.9a) 

0P^(x,t) := h(k-x). (4.9b) 

Note that 

VuPW^ J;2^/,(k-x,t), V0P'^ = (fci,A:2)h,(k-x). 


^We may let h = h(z, t) a priori, but our ansatz for the form of solutions will force dth = 0, so we assume 
that h is time-independant from the outset. 
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This implies V ■ = tr(Vu^^) = {ki + k 2 )fz = 0, i.e., is divergence free. 

Furthermore, 

■ V)u^^ = (1, l)(A;i +/ca)//^ = (0,0) and (4.10) 

(u^^^ ■ = (fci + k2)fK = 0. 

Assume that and 6^"^ defined as in fl4.9p give a solution for some pressure p. 
Substituting these relations into fll.lal) and fll.lbl) . with z = k ■ x and f{xi,X2,t) = 
f{z, t), we obtain 


dtf +Pxi = 

dtf +Px2 = i^\k\‘^fzz +gh. 


(4.11a) 

(4.11b) 


Applying dx 2 to fl4.11al) and dx^ to fl4.11bp and subtracting the results, we obtain 

{ki - k2)dtfz = T^\k\^{,ki - k2)fzzz + hghz, (4.12) 


Since ki + k 2 = 0, it follows that fci — fca = 2fci. Dividing fl4.12p by 2ki and integrating 
in 2 ; implies dtf = + ^gh + ip, for some undetermined fj = Since / 

and h are mean-free and periodic, we conclude that -0 = 0. Thus, / satisfies the heat 
equation 


dtf = J'|kp/« + yii. 


(4.13) 


As usual, we can solve for the pressure by applying dx^ to fl4.11ap and dx 2 to fl4.11bp 
and adding the results to obtain 


Ap(x,t) = gk 2 h^{k ■ x,t). 


(4.14) 


This, along with the mean-free and periodic boundary conditions, defines the pressure 
uniquely. In fact, we can get a more explicit formula for the pressure gradient by 
subtracting fl4.11bp from fl4.11ap to obtain Px 2 — Pxi = gh. Adding fl4.11ap to M.llbp 
and comparing with fl4.13|) yields Px^ + Px 2 = 0. Thus, 


Pxt = -\gh and Px 2 = \gh. (4.15) 

Conversely, let h = h{z) be any mean free, periodic function and let / satisfy 04.131) . 
Dehne (u, 6) by 04.91) and the pressure by 

p = p^^ ■= —^^gH(k ■ x), where H'{z) = h{z). (4.16) 

ZiKi\ 

It is easy to see that (u^^, 6^^) is a solution of (O. We note that the solutions 
of the form (u^^,^^^) thus dehned are distinct from (u'^,6'^) and (u'^,6'^), e.g., 
since both components of are non-zero. Moreover, they are steady states if / is 

15 


































time-independent and is given by / = f steady, where f steady = f steady {z) is the unique, 
mean-free periodic solution of the equation 

J^fsteadyiz) = (4.17) 

Steady State: Due to 04.161) and 04.131) . it is easy to see that (u^^, 0^^) approach a 
steady state which is a plane-wave solution corresponding to the same h and / = f steady 
where f steady is as in 04.171) . 

Remark 4.1. The solution families (u^, 6 '^), and given above have 

unique steady states which they approach exponentially fast in L^. Since these steady 
states are determined by 6 ^, 9^, and h, respectively, which are arbitrary smooth, mean- 
free, periodic functions, this show that the set of steady states of system (Ea is fairly 
rich, and in particular, it is both infinite-dimensional (in the sense that it contains infinite 
dimensional subspaces) and unbounded. 

4-2. Weak Sigma-Attractor 

Definition 4.2. Let S{t) be the semi-group operator associated with 01.ip via 03.14p . We 
dehne the weak sigma-attractor A, of the semi-dissipative system 01 .ip to be the set of all 
(uq, 9q) E H X with the property that 

(i) There exists a global trajectory {u(t),9(t)) dehned for all f G M such that (n(t),6'(t)) 
belongs to H x Lf and solves 01.ip for alH G M and moreover (u(0), 9(0)) = (uq, 6*o). 

(ii) The trajectory (u(t),9(^)) is globally bounded, i.e., the set {(u(t), 6 *(t)) : t G M} is 
bounded in 77 x L^. 

Clearly, it follows from the definition that every bounded, global trajectory {(u(t), 9(t)) : 
t G M} lies entirely on the attractor. The justihcation for referring to A as the weak-sigma 
attractor will be provided towards the end of this section. We will show that .4, is a non¬ 
empty, proper subset of the phase space. 

We will now justify our terminology, namely that of the weak sigma-attractor, for A. 
The global attractor for a dissipative system is often dehned as the cj-limit set of the com¬ 
pact absorbing ball. However, system fll.ip does not have a compact absorbing ball. Yet, 
since || 6 '(f)||z ,2 = HT'oll L 2 is time-invariant, equation fll.lap . considered alone, has a compact 
absorbing ball in 77 by the standard theory of the 2D Navier-Stokes equations (see Ap¬ 
pendix IHU]). By varying the value of IIT'oUlz, one obtains a family of absorbing balls. This 
observation is exploited in Section [5] to write A as the union of cu-limit sets. 

Given initial data (uq, 9o) E H x Lf ^ dehne the dimensionless, time-independent Grashof- 
type number for fll.ip to be 


G := 


5'II^o||l2 




2 

0 


(4.18) 
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Standard energy estimates yield 


||u(t)||i2 < e-""o*||uo||i2 + (4.19) 

Thus, there exists a time t* = t*(||uo||L2) such that, for t > t*, ||u(t)||2,2 is in the ball of 
radius 2i>G in H. For example, t* can be taken as 

^*(||uo||l 2) = ^max|l,log|^^| . (4.20) 

Definition 4.3. Let us denote the following Cartesian product of balls: 

Br := {(uo,6'o) e H X : \\9o\\l^ < r, ||uo||l 2 < R{r)} , R{r) := 2^. 

VKq 

Due to (12.1411 and fl4.19p . for a hxed r > 0, the set is semi-invariant for all positive 
times, i.e.. 


S(t)Br C Br for all f > 0, (4-21) 

where S(t) is the solution semigroup for (11.ip dehned in (13.1411 . 

Definition 4.4. The local attractor at level r, denoted by Ar, is dehned to be the w-limit 
set of Br] that is. 


Ar := u;(Br) := f| '^^{S(t)(uo, 9o) ■ t > r, (uoAo) e SJ, (4.22) 

T>0 


where the closure is taken in the weak topology oi H x L^. Note that by the Banach-Alaoglu 
Theorem, Ar is weakly compact in iL x Furthermore, by Proposition 15.31 (in the next 
section), for each hxed t, the map S{t) : Br —> Br is weakly continuous (i.e., continuous 
with respect to the weak topology on Br). Thus, 

'^\S{t){no,eo)-.t>T, (uo,0o)e54 = S{t)Br. (4.23) 

We will show later that the weak sigma-attractor A (dehned in Dehnition 14.21 as the set 
of trajectories which are uniformly bounded in iL x for all f G M) is indeed equal to the 
union of the sets Ar, i.e.. 


^ = IJ = IJ Ar, (4.24) 

r>0 r>0,rGQ 

where the last equality follows from the fact that the sets Ar are increasing in r. Since Ar 
is weakly compact in if x and .4. is a countable union of weakly compact sets Ar, this 
justihes our referring to A as the weak sigma-attractor. We will show later that in fact A 
(weakly) attracts all bounded sets. 
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Let us recall that the weak topology on a separable Hilbert space is metrizable on 
bounded sets. In fact, one can dehne a metric d which is independent of the bounded 
set so that the corresponding metric space topology coincides with the weak topology on 
every bounded set. This metric can be dehned as follows. Fix a countable dense 

subset of the aforementioned Hilbert space, and let X be a bounded set in it. Dehne the 
metric d on X by 





(0 — ijj, (j)j) I, for all 0, -0 G X. 


(4.25) 


This metric space topology coincides with the weak topology on X. Henceforth, d will 
represent this metric on any bounded set. Due to this fact, namely that on bounded sets, 
the weak topology coincides with the one given by the metric fl4.25p . we have 


dfl) . 3 tn ^ OO and (ug ^ By 

such that S'(f„)(uo,„,6'o,„) (ua,^^)}. (4.26) 


Remark 4.5. In view of the backward uniqueness of sufficiently smooth trajectories, it is 
natural to ask whether there exists a global attractor in a stronger space, such as in H x L^, 
or 1/ X or another space where global well-posedness holds (cf. Theorems 12.31 and 12. 4|) . 
In the case of the 2D Navier-Stokes equations with time independent force in H, it is known 
that the glo bal attractor in the phase spaces V and H coincide and moreover, it is contained 
in [3^. However, unlike in the case of the 2D Navier-Stokes equations, it is not 

known whether one can bound u in uniformly in time. Without such a bound, we are 

not able to prove the existence of a stronger global attractor than the one described here. 


5. Structure of the Attractor 

In this section, we will state and prove the main result describing the properties of the 
attractor. We will need the following dehnition for the semi-distance between two bounded 
sets in the weak topology. For two bounded sets A,B, we dehne the semi-distance 


dist(A, B) 


sup ini d{x^y) 

x&A 


(5.1) 


It is easy to check that for two bounded sets, dist(A, H) = 0 if and only if A G B. 

We are now ready to state our main theorem summarizing the properties of the weak- 
sigma attractor attractor dehned in Subsection 14.21 


Theorem 5.1. The global attractor A defined in Definition \4 .has the following properties. 


(i) The relation A= \J Ar = IJ Ar holds, where Ar is defined by 04.2211 . 

r>0 r>0,rEQ 


(%t) A contains all the steady states of 01 .ip . 
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(in) A has empty interior in the strong (and therefore weak) topology of H x . 

(iv) The weak sigma-attractor A of the system fll.ip is a non-empty, proper subset of the 
phase space X := H x which moreover contains infinite-dimensional sub spaces 

of the phase space. 

(v) The set A is a-compact. More precisely, A is a countable union of weakly compact 
sets Ar, r G Q. 

(vi) The set A attracts all bounded sets. More precisely, if X G H x is a bounded set, 
then there exists r > 0 such that dist{S{t)X, Ar) —)■ 0 as t —)■ oo. 

(vii) For each r > 0, the set Ar, as well as the set A, is weakly connected. 

(via) The global attractor A is invariant (i.e. S{t)A = A) for all t > 0. Moreover, it is the 
minimal set that attracts all bounded sets. 

(ix) (Tracking property) Let {u{t),9{t)),t > 0 be a trajectory in H x There exists 
a global trajectory {uoo{s),9oo{s)), s G M included in A, with the property that given 
e,M > 0, there exists T > 0 satisfying 

sup {||u(s + T) - Uoo(s)||i 2 + ||6'(s + T) - 6 'oo(s)||h-i} < e- (5.2) 

Before proceeding to prove the main theorem, we need to establish the following crucial 

lemma. 

5.1. Auxiliary Lemma 

Lemma 5.2. Let (uo^„, ^ and t^ ^ oo be such that 

S{Q{uo,n,0O,n)^{na,9a) & H X L\ (5.3) 

Then (u^, 9a) G A. More precisely, there exists a complete bounded trajectory (uoo(t), 9^{t)) G 

Br for all t eM., such that 


(Uoo(0),6'oo(0)) = (Ua,^^). (5.4) 

Moreover, G Id and there exists a dimensionless, absolute constant C > 0 such that 

||uoo(^)||hi < V f G M. (5.5) 

Proof. We will start by recalling certain a priori bounds; see [sl and [^. Let (u(f), 9{t)),t > 
0 be any trajectory starting in Br and observe that {u(t),9{t)) G Br due to the semi¬ 
invariance of Br. In particular, this means 

||u(t)||L2 < R{r) and ||6'(f)||Hi <r {t > 0), (5.6) 
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(5.7) 


Moreover, the momentum equation fll.lajl readily yields the bound 

I|u(t)||i 2 + u f ||u(s)||hi ds < ||uo||i 2 + (0 < t < T), 

Jo 

where the invariance of the norm of 6{-) along a trajectory is also used. Thus, there 
exists t G (0, T) such that 


u 


Hi — 


< 2 


|uo|li2 ^ 


oiii. 


z/T 


2^2 




(5.8) 


Since the distribution function of 6 remains invariant along a trajectory forward in time, by 
replacing the initial data if necessary, we can (and henceforth will) without loss of generality 
assume that uq G Id and moreover, setting T = in flS.Sp . we may also assume 


uoIIhi < 2i?(r)Ko, 


(5.9) 


where R{r) is as in Dehnition fl4.3p . With the assumption fl5.9p . noting that the periodic 
boundary condition and fll.lcp implies (in 2d) that (i?(u, u), ^du) = 0, and consequently, 
following techniques in (Chapter XIII, pages 111-113), we additionally have the uniform 
bound 


sup ||u(t)|lHi < CR{r)Ko, ( 5 . 10 ) 

t>ti 

where R{r) is as in Dehnition 14.31 ti = Additionally, from the hrst relation in flh.lOp 
and using again (i?(u, u), Au) = 0, for any r > 0 and T > 0, we have 

l|u(r+ T)||^i+2z/^ ||Au||^2 < C^R(r)^K^+ ^ 11^011^2 ^^ 

Let now (uo,n, 6'o,n) £ Br and the sequence of times ^ oo be such that fl5.3p holds. 
Translating in time by G [0, ,if necessary, we will assume without loss of generality that 

uo,n £ y, and that the bound in fl5.9p holds for all uo,„. Since uo,n W, h immediately 
follows that 

UaEV and ||ua||ei < 2i?(r)Ko- (5.12) 

Set 

(u„(t), Onit)) := S{t + t„)(uo,„, 6o,n), t G [-R, oo). (5.13) 

Note that the trajectories (u„(t),6*n(t)) are dehned, for n sufficiently large, on the intervals 

f M M\ 

Im-=( -, MgN. ( 5 . 14 ) 
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By fl5.6jl and fld.lOp . for n sufficiently large, the pair (u„, 6 *„) is bounded uniformly in 
V)xL°°{Im, L'^), and moreover due to fid.lip . u„ is bounded uniformly in EI^)- 

The above mentioned bounds are uniform with respect to n, although they may depend on 
M in general. From the functional equation fl2.16p . the time independent bound fl5.6p 
on 9n{t), the (uniform in n) bound on in H^), and the Ladyzhenskaya inequality 

02.41) ■ we readily obtain that is uniformly bounded in Consequently, {u„} 

is an equicontinuous family (in time) in Moreover, from the Agmon inequality 

02.51) . and the uniform bound on u„ and Ol.lbp . it also follows that ^9n is uniformly 
bounded in (in fact, one can show that the bound is uniform in 

for any p < oo). Thus, {9n} form an equicontinuous family in Now, by 

the Arzela-Ascoli Theorem and the Cantor diagonalization procedure, there exists a pair 
(uoo(-)) ^oo(-)) ^ C{{— 00 , 00 )] H) X C((—cx), cxo); and a subsequence (which we also 
denote by (u„, 9n)) such that for all M E N, 

(u„(-), 0n(-)) ™ (Uoo(-), ^oo(-)) in C{Im] H) X C{Im] H"'). (5.15) 

Note that since Un(0) is equal to the u component of S'(fn)(uo,n, 6 ^ 0 , 71 ), by the limits assured 
in fl5.3p and fl5.15p . we have Uoo(O) = u^. On the other hand, due to the hypothesis and 

9{tn) = 6 * 71 ( 0 ) 9a and 6 *(t 7 i) = 6 * 77 ( 0 ) 6 * 00 ( 0 ) as n -)■ 00 . (5.16) 

Thus, 6 * 00 ( 0 ) = 9a. 

Next, we will show that, by passing through a further subsequence if necessary, we can 
ensure that for all rj G the functions ( 6 *oo(-))h) ^ C'((~oO) 00 );®)) nnd that 

( 6 * 77 (-),h) ^ ( 6 *oo(-)>h) in C(/m;M), as ncx), (5.17) 

and also that (uoo(t), 6 * 00 (t)) E Br'i t E M, i.e., 

2 2 

|| 6 *oo( 6 )||i 2 < and ||uoo(t)||i 2 < 4^^ V f G M. (5.18) 

To see this, note that for any G (which, in 2D,implies that VC G -h°°), we have from 

dnbi), 

I ( 6 * 77 ( 0 ) - 6 ', 7 (fi),C) I < [ I ( 6 'n(s)u„(s), VC) Ms (5.19) 

Jti 

< ||6*77(-)lln°°(/M;i^hll^"(')lln°°TM;^^)II^CIU°°|0 - Ol- 

Thus, for all C G the functions {( 6 *, 7 (-), C)} equicontinuous on Im- Due to the uniform 
bound on || 6 * 77 (t)|M 2 , we have a pre-compact family. Let {Cj}^i be a set of functions in 
which is dense in . By the Arzela-Ascoli Theorem and the Cantor diagonalization process, 
we may also assume that 

(^n(-),0) ^ (^oo(-),0) in C(Jm;M) for all M,j E M. (5.20) 
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To complete the proof, we show that {9n{-), (} —>■ (9oo(-), C) uniformly on Im for any ( G L^. 
Let e > 0 and choose j such that |C ~ 01 ^ 2 p^jrT- Then, 


I (9n(t) - 9m(t),0 I < I (9n(t) - 9m(t),Q} I + Pnit) - 0,„(t)|UOIO " CIU^ (5.21) 

< I {9n{t) - Cj) I + 2||6'o||l2 -. 

2||c'o||l2 

This shows that {( 6 '„(-), C)} is uniformly Cauchy in C{Im)- Since 9n{-) converges to 9oo{-) in 
the uniform limit of {9n(t)X) is {9oo(t),C)- Due to weak convergence of 9n(t) 
and strong convergence of u„(f), the uniform bounds fl5.18p readily follow. This hnishes the 
proof of the claim. 

Next, we will show that (uqo, 9oo) is a global solution of fll.ip . i.e., we need to show that 
(uoc^'oo) satishes fl2.9p for appropriate test functions $ and 0 (see Dehnition 12. ip . Note 
that {Un,9n) satisfy the weak formulation fl2.9p . Passing to the limits in the linear terms, 
due to fl5.15p . we hnd 


(u„(s),<h'(s)) ds ^ / (Uoo(s),<h(s)) ds. 


u / ((u„(s), <h(s))) ds^u ((uoo(s), <h(s))) ds, 


{9n{s)g,^{s)) ds -)■ / ( 6 'oo(s)g, $(s)) ds, 

do 

{9n{s),(j)\s)) ds ^ [ {9oXs),(j)'(,s)) ds. 


It remains to show the convergence of the remaining non-linear terms. Let 


J(u) := / (u„ ® u„ - Uoo ® Uoo, v<h(s)) ds. 


J{n) := / (u„(s) 6 'n(s) - Uoo(s) 6 'oo(s), V0(s)) ds. 


(5.22a) 

(5.22b) 

(5.22c) 

(5.22d) 

(5.23a) 

(5.23b) 


Jo 

The convergence I{n) —)■ 0 as n —)■ cxd follows from the convergence of the component 
guaranteed in fl5.15p . To show J(n) —)■ 0 as n ^ oo, we write J{n) = Ji{n) + J 2 {n), the 
dehnitions of which are given below. We have 


Ji{n) := f ((u^(s) - Uoo(s)) 6 'n(s), V0(s)) ds-)■ 0 (5.24) 

do 

as n ^ oo, since u„ —)■ u in C{Im, H) and is uniformly bounded in H). For J 2 , 

we have 

J 2 in) := f (uoo(s)(6'n(s) - ^ 00 ( 5 )), V0(s)) ds-)■ 0, (5.25) 

Jo 
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due to the fact that 0 is a smooth test function for each hxed M, sup^g/^^^ ||uoo(s)||hi < oo 
and 9n{-) 0oo{-) in 

Since (uoo(t), 6 *oo(t)) is a global trajectory with 

(uoo(t), 6 *oo(t)) G-Br for alH G M, (5.26) 

it readily follows from flh.lOp that ||Uoo(t)||Hi < CR{r)K,o for all f G M. □ 

We will also need the next proposition in order to proceed. 

Proposition 5.3. For each fixed t and a sequence {uQ^n,9o,n) (uo,6'o), we have 

5(t)(uo,„,V) ^>5(t)(uo,0o). (5.27) 

In particular, for each fixed t, the map S{t) : X —> H x is weakly continuous, where X 
is a bounded subset of H x . 

Proof. Fix T > 0 and let (uo,n, 9o,n) be the sequence in the hypothesis of the proposition. 
Due to the weak convergence, the sequence is uniformly bounded in the phase space. Let 
fvinit), 6n{t)) := S'(t)(uo,n, 6 *o,n), 0 < t < T. By arguments similar to the proof of Lemma 
15.21 and by the uniform (in n) bounds provided in Theorem 12.21 (u„(-), 6 *n(-)) converges 
(in C[e,T]]H x H”^) for all e > 0) to a solution (u(-),0(-)) with initial data (uo,0o)- By 

uniqueness of solutions, {u{t),6{t)) = S{t){uo,6o),0 < t < T. Since (u„(T), ^^(T)) 
(u(T), 9{T)) (see proof of Lemma [5.2p . the proof is complete. The weak continuity of S{T) 
on bounded subsets now follows from noting that the weak topology is metrizable on bounded 
sets. □ 

5.2. Proof of Theorem \5.1\ 

For the remainder of this section, we focus on proving Theorem 15. II In particular, we will 
examine to what extent A shares the properties (jl]), (ini), and fillip with .AnS) the attractor 
for the 2D Navier-Stokes equations. 

Proof of Theorem 15.11 (i) 

Proof. We first show the inclusion in the “d” direction. Let (u^o, doa) G Ar for some 
r > 0 . Then, since Ar = oj{Br), due to fl4.26p and Lemma E21 it immediately follows that 
(H-00 5 ^Oo) G a. 

Next, we show the inclusion in the “c” direction. Choose any (uq, Oq) G A, and suppose 
that (uo, ^o) ^ As for all s > 0. By the dehnition of A, there exists a trajectory (u(t), 9(f)) G 
A such that for some t^, (^u(to),9(fo)) = (uo, 6 *o). Furthermore, (^u(f),9(f)) is bounded in 
H xL^ uniformly for all f G M, so there exists an r > 0 such that (u(t), 9{t)) G Br for all f G M. 
Since (^u(f),9(f)) is dehned for f G M, we must have (uq, 9o) = S(T)(u(to — r), 9(to — T)) G Br 
for every r G M. But, (u(to — t), 6 *(to — n)) G Br, so (uq, 6 *o) G u}{Br) = Ar- □ 

Proof of Theorem 15.11 (ii) 
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Proof. Since any steady state is bounded for all times t G M, the proposition follows imme¬ 
diately from the definition of A. □ 

Proof of Theorem 15.11 (hi) 


Proof. If (Ua, 6 a) G A, then by Lemma [5.21 G V. Any open ball in if x (in the norm 
topology) must contain points whose hrst component is not in V. Thus, there does not exist 
an open ball contained in A. □ 

Proof of Theorem 15.11 (iv) 

Proof. Theorem 15. II fiiil immediately implies that A^ X. To show it is nonempty, consider 
a solution of the form ,U 2 ,6^,p^), defined in Subsection 14.11 Let us also impose that 
dtU 2 = 0, so that = U 2 , where U 2 = U 2 {xi) is the unique, mean-free and periodic solution 
of the equation 

(P 

h^M2(xi) = -g9^{xi). (5.28) 

Then (wf^, U 2 , 6^,p^) is a steady state solution of fll.ip . so {uf, U 2 , 6^) G A. Since 9^ G 
can be chosen arbitrarily, A contains an infinite dimensional subspace. □ 

Proof of Theorem 15.11 (v) 

Proof. Since A^ C Ar 2 if < r 2 , it follows that A = IJ Ar- Po complete the proof, it 

r>0,rEQ 

is enough to show that for every r > 0, is weakly compact in ii x Lf. By the Banach- 
Alaoglu Theorem, 5^ is weakly compact in ii x Lf. Since S{t) is weakly continuous, S{t)Br 
is weakly compact. Thus, Ar := u{Br) is the intersection of weakly compact sets, so it is 
weakly compact. □ 

Proof of Theorem 15.11 (vi) 

Proof. Suppose X C ii x is bounded. Then there exists r > 0 such that X G Br, where 
Br is the absorbing set as defined in Definition 14.31 Thus, due to semi-invariance of Bj. under 
S'(i), we have S(t)X C B^. for all t > 0. The proof now proceeds by contradiction. Assume 
that the conclusion of the theorem is false. Then there exists a sequence (u„, On) e X, 
tn ^ 00 and e > 0 such that 

d{S{tn)iUn,9n),y) > 6 for all y G u{Br) = Ar G A. (5.29) 

However, since S'(i„)(u„, 9n) G Br, by the weak compactness of Br, there exists yo G Br and a 
subsequence Uj such that ^(i^J(u„., yo, or equivalently, d{S{tnj){Un.,9nj), yo) 0 . 

By the definition of oj{Br), the point yo € oj{Br) and this contradicts fl5.29p . □ 

Proof of Theorem 15.11 (vii) 
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Proof. Suppose Ar is not weakly connected. Then there exist open sets Oi and O 2 in the 
weak topology oi HxL"^ such that Ar C O 1 UO 2 , 0 in 02 = 0, ArPOi 7 ^ 0, and ArP[02 7 ^ 0- 
Recall that Ar := uj{Br). Since Br is weakly connected and S{t) is weakly continuous, S(t)Br 
is also weakly connected. Thus, for each n G N, there exists (u„, 6*„) G {S{n)Br) \ (Oi U02)- 
Since Br is invariant under S(t), (u„,6*„) is bounded. Thus, by the Banach-Alaoglu 
Theorem, there exists a point (u, 6*) ^ H x and a subsequence ^ (u, 6*) in the 

weak topology of H x Since Oi U O 2 is weakly open, (u, 9) ^ OiU O 2 . But, by Theorem 
15.11 Ivih Ar attracts all bounded sets in the weak topology, so that the weak limit of (u„, 9n) 
must he in C Oi U O 2 ; a contradiction. Therefore, Ar is connected. 

To show that A is weakly connected, note that by Theorem I5.1f i). A = [Jr>QAr. Fur¬ 
thermore, each Ar is connected and contains the zero element. Therefore, if Oi and O 2 are 
weakly open sets that separate A, all Ar must be contained in the one containing the zero 
element, so A must be contained entirely in either Oi or O 2 , so that A is connected. □ 

Proof of Theorem 15.11 (viii) 

Proof. Let (uo,6'o) G A. Then there exists a bounded, global trajectory {u{t),9{t)),t G 
M passing through (uo,6'o). It is clear from the dehnition of A that all points on this 
trajectory also he on A. In fact, since this trajectory is bounded, it lies entirely in some Br 
and consequently, in u{Br) = Ar- We may assume, by translating time if necessary, that 
(u(0),6*(0)) = (uo,6*o). Then, by (forward in time) uniqueness of solutions, 

9{-t)) = (u(0), ^(O)) = (uo, 6»o). (5.30) 

As remarked above, (u(—t), 6*(—t)) G A. Thus, A C S(t)A for all t > 0. For the reverse 
inclusion, simply note that the (forward in time) uniqueness of solutions guarantees that 
{u{t),6{t)) = S{t){uo,6o) for all t > 0, where {u{f),d{t)),t G M is the above mentioned 
bounded global trajectory. Thus S(t)A = A. 

The fact that A attracts all bounded sets was proven in Theorem 15.11 (vi). We will now 
prove that it is minimal. Let A! be a set which has the property that it attracts all bounded 
sets and let (uq, 6q) G A. Then there exists a bounded, global trajectory passing through it. 
In view of the previous paragraph, we see that the set in the phase space 

g = {(u(t), e{t)) : t G M} C FT X (5.31) 

is invariant under the semigroup S{t) for all f > 0, i.e., S{t)g = Q for all t > 0. On the 
other hand, ^ is a bounded set in iF x F^. Thus, by the assumption that A! attracts all 
bounded sets, 

limsup dist(S'(t)^, ^') = 0. (5.32) 

t—yoo 

□ 


Thus, Q C A', and consequently, (uq, 6q) G A!. 
Proof of Theorem I5.1f ixj 
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Proof. Note first that since (u(t), 9(t)),t > 0 is bounded in H x Lf, there exists (ua, 6a) and 
a sequence such that (u(t„), (uaj^a). By Lemma E21 we immediately infer 

that (Ua,6'a) € A. Furthermore, the global trajectory (uoo (-)5 ^oo(-)) constructed in Lemma 
15.21 lies in A and 


{\ln{t)An{t)) ■= + + G [-tn,Oo), (5.33) 

converges to (Uoo(-)) ^oo(')) x C([—M, M]; EI"^) for all M > 0. This 

concludes the proof. □ 


6. The Presence of 2D Turbulence 


The Batchelor-Kraichnan-Leith theory 0, 35| of 2D turbulence (inspired by that of Kol¬ 
mogorov in 3D jH, HH, 11^) asserts that, on average, the behavior of eddies in turbulent flows 
is determined by their length scales. In a relatively large range of scales [k*, k*] (called the 
inertial range) viscous effects are negligible and enstrophy is transferred at a nearly constant 
rate rj from one length scale to the next smaller one (termed the enstrophy cascade). The 
dissipation range consists of the very small length scales where the viscosity annihilates the 
enstrophy. Heuristic arguments by Batchelor and Krachnan j^, [s^ place the dissipation 

1 If? 

range beyond a wave number where p is the average rate of enstrophy dissi¬ 

pation per unit mass. The main tenets of this empirical theory can thus be summarized as 


(see, e.g. 


0 , 00 , 0 ), 


a signihcant amount of enstrophy is in the inertial range [k*, k* 


(ii) this range is wide, i.e. k* << k* ~ Krj] 

(hi) the direct cascade of enstrophy (to smaller scales) holds over this range; 

2/3 

(iv) the power law eK, 2 K ~ holds for the amount of energy eK, 2 K contained in the length 
scales K to 2k, for k E [k*, k*]. 


Rigorous justihcation for parts of the theory has been obtained in the series of works 


121, ll3, l2l|, [2^ among others. In |l6l.ll8l|. it is shown that many of the ubiquitous averages 


in empirical turbulence theory, upon the application of which patterns are observable, can 
be taken to be hnite time averages, albeit on sufficiently long periods of time. Furthermore 
in 0 extension of the above mentioned works has been obtained for forcing in all scales. 
Thus, in view of the results in 0 16, 181, the universal features of turbulence already hold 
for hnite time averages of the form 

(d>) := 


ri2 


to -t 


1 Jtl 


<F(u(s))(is, 


( 6 . 1 ) 


where $’s are the relevant physical functionals on the phase space H and ^2 > H > 0 satishes 

G 


max{ti, t 2 —ti} » 


UK, 


2 • 
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Here G is the Grashoff number which is a non-dimensionalized version of the norm of 
the driving force in the Navier-Stokes equations. 

In principle, the results in apply in our setting to the velocity equation written in 
the functional form in fl2.16p . However, for the inertial range to be sufficiently large for 
turbulent patterns to emerge, it is necessary (though not sufficient; see Remark 16.21 below! 
for the magnitude of the driving force, as measured by the Grashof number, to be large. 
Indeed, it is well-known that the 2D Navier-Stokes equations converge to a steady state 
if the force is time independent and the Grashof number is sufficiently small 
Accordingly, let us dehne the dimensionless (time-dependent) number G^ by 


41, 40, 42 


G. = G^{t) : = 


\Pa{g6){t)\\L^ g 


2^2 




0 


2^2 




11^0 


|2 

Il2 


- I|ft29(i)lli. 


( 6 . 2 ) 


The a here is used to denote the fact that the norm of the solenoidal projection of the force is 
taken, rather than the norm of the force itself. Note that, as in 04.181) .one may also consider 
the somewhat simpler (though potentially larger) time-invariant, dimensionless number 


G : = 


fi'll^o||L2 




2 ’ 
0 


(6.3) 


which obviously satishes 0 < Go-(t) < G. The effective Grashof number governing the 
dynamics of 02.161) is dehned to be 


G* = limsup Go-(t). (6.4) 

t—)-oo 

Thus, the complexity of the flow, at least in regard to the statistical features, is expected 
to be determined by G*. This is borne out by the following result which is analogous to the 
2D Navier-Stokes equations. Before stating this result, let us observe that since 9 G the 
distributional derivatives dxi9,i = 1,2 belong to and moreover, 

ll^xi^(^)||H-i ^ II^(^)I|l 2 = ||0o||L2,i = 1,2. (6.5) 

Theorem 6.1. Assume that G* < e where e > 0 is sufficiently small and that P„g9(t) 
converges weakly to f G H as t ^ oo. Then u(t) approaches a steady state of the 2D Navier- 
Stokes equations with time-independent driving force f. Consequently, if dxj^9(t) converges 
to zero in as t ^ oo, or equivalently, if \\Pfj{g 6 ){t)\\L 2 converges to zero, then u(t) 
converges to zero in H. 

Proof. For notational simplicity, denote F(t) = P^^gOf) and observe that F(t) G FT and 
||F(i^)||L 2 < g\\ 0 o\\L 2 . Moreover, since F(t) ^ f weakly, we also have 

||f||i 2 < limsup ||F(t)|| 2,2 < (6.6) 

>-oo 

Thus, by shifting time if necessary, we may assume without loss of generality, that 

max{||f||i 2 ,sup ||F(f)||i 2 } < 2ez/\o. (6.7) 

t>o 
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Furthermore, due to 05.101) in the proof of Lemma 15^ and by shifting time again if necessary, 
we also have the uniform bound 


sup ||u(f)||Hi < Cez/Ko, (6.8) 

t>o 

where C is a non-dimensional, absolute constant. Recall also that on bounded subsets of 
H, the weak topology is metrizable and a bounded sequence v„ —)■ v weakly if and only if 
- v)||x ,2 -)■ 0. Thus, 


lim _f)|| _ Q_ 

t—>-oo 


Let v(f) be the solution of the Navier-Stokes equation with initial data uq = 0, and the 
force given by f. Denote w = u — v. Then w solves 


dw 

dt 


uAw -|- B{u, w) -|- B(w, u) = F(t) — f, V ■ w = 0, w(0) = uq. 


( 6 . 10 ) 


Taking inner product of the above equation with w and by applying Young’s inequality, we 
readily obtain 

^^I|w||i 2 + ^\\A^w\\l 2 < ^\\A-^F{t) - i)\\l 2 + I (R(w,u),w) |. (6.11) 

From fl2.8b[) . we get | (R(w,u),w) | < ||y4^u||i2||-w||^2||RW||j;,2. Thus by Young’s inequality, 
we have 


KF(t) - f)||i2 -h -||^2 u||^2||w||^2 

< ipn(F(«) - f)ii|, + c■v..K?llw||^ 


( 6 . 12 ) 


where the inequality in the second line is obtained using 
inequality, we obtain 


If CY" < i by Poincare 


— ||’w||^2 ^^||w||^2 < -\\A 2(F(f) - f)||i2. 


(6.13) 


By Gronwall inequality and fl6.9p . it immediately follows that limj_!.oo ||w(f)|| = 0. On the 
other hand, it is well-known that if ||f||L 2 < ez/^^Q and e is sufficiently small, then the solution 
to the Navier-Stokes equations corresponding to the time-independent force f, converges to 
a hxed point 40|, 1^, |^. Thus, if e is sufficiently small, v converges to a hxed point vq G if 
which satishes uAvq + B{vo, vq) = f. Since w(t) converges to zero in H, it follows that u(t) 
converges to vq in H. 

For the second part of the theorem, observe that by fl2.19p . 


ip.sdmL^ = WRiomL^ = 


(6.14) 


It follows that if dx^9{t) converges to zero in as f —>■ cx3 if and only if ||Po-g0(f)||j;^2 
converges to zero. Thus one can apply the hrst part of the theorem with e = 0 to conclude 
that u(t) —)■ 0 in FT as f —?• oo. □ 
















Remark 6.2. It is useful to examine Theorem 16.11 in light of the examples presented in 
Subsection O In the “vertical solutions” presented in (i), due to fl2.19jl . P^g9(t) = 0 and 
in conformity to theorem 16.11 u{t) —?• 0. On the other hand, for the “horizontal solutions” 
presented in Subsection 14.11 fii). due again to fl2.19p . \\Po-g9(t)\\L2 = || 6 * o || l 2 - Consequently, 
G* can be made arbitrarily large, yet the flow is not turbulent as it converges to a (laminar) 
steady state. Incidentally, in addition to the example provided by Marchioro [^, ^ with 
forcing in the hrst eigenmode of the Stokes operator, this provides another example of 
solution to the 2D Navier-Stokes equations with arbitrarily large Grashof number for which 
the solution converges to a steady state. 


Remark 6.3. If one considers equations fll.lal) and fll.lcl) . this is the usual 2D Navier-Stokes 
system with forcing g9, and so the framework developed recently in for 2D turbulence 
with forcing at all scales applies. Moreover, unlike the usual 2D Navier-Stokes equations, the 
forcing here is not ad hoc, but rather intrinsic to the system. Noting that ||g 6 *|| 2,2 = ||g 6 *o||L 2 
for all f > 0, it seemingly allows arbitrarily large Grashof number as well. However, due 
to the divergence free condition on the velocity held, the effective driving force is given by 
Pfjg9(t). The quantity \\Pa-g9(t)\\L2, although bounded above by ||g 6 *o||L 2 , can potentially 
approach zero for large time. Thus, if one can show that = liminfi_j.oo \\Pag9{t)\\L2 can 
be made arbitrarily large, then one can assert that fll.lap and fll.lcp . or equivalently fl2.16p 
contains the entirety of Grashof numbers involved in 2D turbulent dynamics. Whether or not 
arbitrarily large and a truly turbulent how can be achieved, remains an open question 
at this point. As observed in Remark 16.21 one can indeed achieve arbitrarily large 
the corresponding how however approaches a steady state and is therefore not turbulent. 
The velocity prohle is also laminar; thus it does not yield a turbulent Lagrangian dynamics 
either. 


7. The Energy-Enstrophy Plane 


In this section, we investigate the relationship between the ener gy a nd the enstrophy in the 
conservative set. Here, we follow many of the ideas of 0 , [ill. Ilil. [isj which investigated 
the relationship between energy, enstrophy, and palenstrophy in the attractor of the 2D 
Navier-Stokes equations. 

For (u, 9) G A, let us dehne 


x{t) ■= 


u(^)IIhi 


and 


m ■■= 


u(t)iii2' 


(7.1) 


After some computation, which was carried out in the context of the 2D Navier-Stokes 
equations in j^, it can be shown that 


dt 


2||u||i2 y A 
2||u||i2 Y 1/2 



(7.2) 
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where w := (^A— ^ ^ u. Let us consider the quantities ||u||i 2 and ||u||jji as variables. 

Level sets of x correspond to curves along which ||u||j ^2 and ||u||]hii are parabolically related. 
Due to fl7.2p . if x(t) > g\\9o\\L'2/i^, then x must decrease. Furthermore, due to the Poincare 
inequality, we must always have X(t) > Ai. Therefore, trajectories flow towards the set 

A := |(u, 6 ') : ||u||^i < ||u||^i > Ai||u||^ 2 , || 6 '||i 2 = || 6 'o||l 2 | . (7.3) 

It follows that Al C A. See Figure [T] for a depiction of the projection of A onto the energy- 
enstrophy plane. 



Figure 1: The projection of the attractor lies inside the shaded region. 


Let us consider the upper-boundary of this set; that is, the parabola ||u||^i = ||u|| 2 , 2 . 

Suppose there is a point of (u, 6) which lies on the parabola, and moreover that u = uq 
is independent of time (a priori, 9 may be balanced by a time-dependent pressure). Then, 
since u(t) = uq, we have X{t) = A(0). On the parabola, equation fl7.2p yields 


0 = 2Au - 


u 


Hi 


U 


s9 

u-= 2Au — A(0)u — 


L2 


(7.4) 


SO that, in fact, 9{t) = 9 q. Taking the inner product with u gives 


2||u||ei — A(0 )||u||^ 2 -h 


(g^,u) 


u 


,A(o) = MA 


u 


L 2 


Owing to fl7.5ll . and the fact the u is on the parabola, it follows that 


(7.5) 


u) 


z/ u 


= A(0) = 


u 


Hi 


9\\0\\ 


L 2 


L2 


U 


L2 




(7.6) 


Thus, 


(g0,u) = \\g9\\L2\\u\\L2. 
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(7.7) 



















By the Cauchy-Schwarz Theorem, we must have that 
where c = Applying this to fl7.4p . we have 


is a scalar multiple of u, say gO = cu 


11^2 


224u = A(0)u + 


gO 



u|Ihi I {sO,u) \ 

u|li2 ||u|| 2 , J 


u 


(7.8) 


so that u is an eigenfunction of A with eigenvalue | ( I—+ irTi^l • 

Going back to the equation {0,g9)^ = g9 = cu = c{ui,U 2 )'^, we observe that we must 
have Ml = 0. The divergence-free condition then implies that 8 x 2^2 = 0, so that u ■ Vu = 0. 
Moreover, since g9 = CU 2 , we have 8 x 29 = 0, so that 8t9 = —u -9 = 0. Thus, 9 depends only 
on xi, and again since g9 = cu 2 , so does M 2 - The entire system fll.ip therefore reduces to 
the following relation: 


8 x 2 P - ^ 8 l^U 2 = g9. (7.9) 

Applying 8 x 2 yields 8 ‘l^p = 0, so that the periodic boundary conditions imply that p is a 
constant. The relation now becomes 


—iy 8 lu 2 = g9 = cu 2 . 


(7.10) 


The periodic boundary conditions further constrain c to be of the form 

for some n G N. 


n^TT^ 


c = c„, := u- 


L 2 


(7.11) 


Summarizing these observations, we have the following proposition, which essentially ap¬ 
peared in 10(1 in the context of the Navier-Stokes equations. 

Proposition 7.1. Suppose {u,9) is a smooth, steady state solution of fll.ip which lies on 
the parabola ||u||^i = ||n ||^2 Then U 2 and 9 depend only on xi, and 


U\ — 0 , U 2 Cn9, 


Cn • — 


gL^ 

Z/n^TT^ 


for some n G N 


(7.12) 


and U 2 = U 2 {xi) is an eigenfunction of the operator — 8 \^ (respecting the periodic boundary 
conditions and mean-free condition). 

Remark 7.2. The above proposition implies that steady states on the upper boundary of 
A are discrete points corresponding to geophysical steady states of stationary columnular 
flows. Somewhat more generally, it can be shown that if (u, 9) is a smooth flow for which 9 
depends upon xi alone, and 9x-^ 7 ^ 0 a.e., then (u, 9) is one of the horizontal flows described 
by fl4.5p and 04.61) . 
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8. Open Questions 

In this paper, we have introduced a notion of an attractor for a semi-dissipative system, 
i.e., a system which is a hybrid of parabolic and hyperbolic equations. While in some 
respects, it shares certain properties with the attractor of a dissipative system, there are 
several open questions concerning the attractor A. 

(i) We showed that the attractor A is in some respects, a thin set, i.e., it has empty 
interior. This proof uses the fact that the velocity component of any point on A is 
regular (belongs to V). On the other hand, we saw in subsection 14.11 that A contains 
infinite dimensional subspaces of steady states. In particular, the projection of the 
attractor on the temperature component contains all functions of the type 6 = 6{xi) 
or 6 = 6 {x 2 ) or functions of the form ^(x) = h(k ■ x), for an arbitrary function of one 
variable h. The question of how rich this set can be remains open. In particular, if 
one takes the projection of A onto the temperature component (say P 6 »(u, 0 ) := 9), the 
following questions remain open 

(a) Does PgA have an empty interior in 

(b) Is P 0 A a proper subset of 

(c) Is PqA dense in 

It is curious to note that in all the examples of steady states obtained here, the tem¬ 
perature 0 is a function of only one variable. 

(ii) We showed that although A is unbounded and inhnite dimensional, it is a-compact 

in the weak topology. In particular, it can be written as a countable union of weakly 
compact omega limit sets Ar = uj{Br) where are absorbing, invariant balls as dehned 
in fl4.3p . The question is whether the attractor has a pancake-like structure. In other 
words, is it true that A = a;(i?r), where 

Br := I(uo, 6 'o) : || 6 'o||l 2 = r, ||uo||hi < R{r) = 2 ^ j? ( 8 . 1 ) 

If this is the case, the u-component of u{Br) can be regarded as an attractor at level 
r for the 2D NSE with time varying force 6 whose L^-norm remain hxed at r. The 
temperature component 6 on the attractor is then a rearranged version of the initial 
temperature 6 q. The problem with this picture is that an element of the omega limit 
set is obtained as a weak limit of points S{t){uo^n,do,n) with (uo^^, 6 'o,n) ^ Br- In the 
weak limit, the norms might decrease, i.e., it is possible that 

lim ||S'(f)(Uo,n, 6 'o,n)||L 2 > II lim 5'(f)(uo,n, 6 'o,n)||L 2 - ( 8 . 2 ) 

n^oo n—>-oo 

This might destroy the simplistic pancake structure described above. Whether or not 
this happens, and if it does, what its implication is for the asymptotic dynamics, are 
questions that remain open. Furthermore, although the whole attractor is inhnite 
dimensional, is there any kind of hnite dimensionality in its constituent pieces AR 
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9. Appendix 

9.1. Absorbing Ball 

For the sake of completeness, we include here a proof of the existence of an absorbing 
ball. By taking the inner product of fll.lal) . we obtain the energy estimate 


^^I|u||i 2 + z/||u||2i = (P^(0g),u) < ||P,(g0)||i2||u| 


L2 






L2- 


(9.1) 


Using the Poincare inequality fl2.ip . fl2.19p . and fl2.14p . we obtain 

d llD^„nMl2 ||fl_||2 


so that 


„I|2 , .„_2||„||2 , WP.ismh Mh-WMh 


'-0 


(9.2) 


u(«)|li= < e-“S'||u„||i, + (||9„||i, - IlfljOolli,) (1 - e-<‘). 




(9.3) 


Thus, 


limsup IIu 

t—)-oo 


V Kq 


(9.4) 


In particular, there exists a time f* = f*(||uo|| 2 , 2 ) such that, for t > f*, u(f) G B 2 ug%, the 
ball of radius 2i/G* in H. For example, f* can be taken as 


u(I|uo||l2) 


1 

—K max 


l,log 


ll^olli2 ] 

3z/2G*2 J • 


(9.5) 


9.2. Mean Preservation 

We show that the mean is preserved by the flow. Since the equation is globally well-posed 
and possesses higher-order regularity, we may assume our solutions are smooth enough for 
the operations below to be justihed rigorously. We give the argument only in two dimensions, 
but it extends without difficulty to the higher dimensional case. 

Integrating fll.lbl) over kl, we hnd 


d 

dt 



V ■ (u^) dx = 0 


due to the periodic boundary conditions. Thus, 


(9.6) 



(9.7) 


33 



















for some dimensional constant ci. Integrating fll.lal) over we find 

d 


dt 


U(ia;+ / V ■ (u ® u) dx + / Wpdx= / z/Auda;+ / gOdx, 


in 


so that 


[ udx = 

dt Jr, 


' dx, 


Writing this in component form, this becomes 

jL 

dt 

jL 

dt 

Thus, 


Ml dx = 0, 


U 2 dx = g / 6 dx = gci- 


U 2 dx = gcit + Co¬ 


rn 

(9.9) 

(9.10a) 

(9.10b) 

(9.11) 


Taking the inner product of fll.lap with u, and of fll.lbp . using fl2.7p . and adding the 
results, we obtain 


liMlh + '^IIuIIhi = (g^,u) < c/||d||i2|u||i2 (9.12) 

Applying Poincare’s inequality and the Cauchy-Schwarz and Young’s inequalities, we find 


1 d 

2 dt II ‘^IIl^ 


u||f 2 + Z/Al||u|| r 2 < -7^ 




+ 




u 


L2- 


Thus, 

|||u||i2 + ^Ai||u||i2<£||d||i2 = pi^||do||i.. 

Integrating this equation, we have 

pt 


|u(t )||22 Y e ^^^*||uo ||^2 + 




_ -J/Aii|| ||2 I 


^0 
1 — e 


2u\\ 


II^o||l 2 ds 


(9.13) 

(9.14) 

(9.15) 


—uXit 


V 


2A? 




o|Il2- 


In particular, u G L°°{0,T; H) for all T > 0. This implies 


l^fCit + Col = 


U 2 dx 


< / |M2Ma; < |r2|^'^^||u||i;^2 


(9.16) 


, , 1 — 


2i/2a; 


0IIl2 


Sending t —)■ cx), we see that we must have ci = 0. Thus, fl9.10ap and fl9.10bp become 


jL 

dt 


udx = 0, 


jL 

dt 


6 dx = 0. 


(9.17) 


Therefore, the flow preserves the spatial mean. 
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